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The traveling salesman problem for lines and rays in the plane 11 
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£SJ ' Abstract 



In the Euclidean TSP with neighborhoods (TSPN), we are given a collection of n regions 
(neighborhoods) and we seek a shortest tour that visits each region. In the path variant, we 
seek a shortest path that visits each region. We present several linear-time approximation 
algorithms with improved ratios for these problems for two cases of neighborhoods that are 
(infinite) lines, and respectively, (half-infinite) rays. Along the way we derive a tight bound on 
the minimum perimeter of a rectangle enclosing an open curve of length L. 
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1 Introduction 



In the Euclidean Traveling Salesman Problem (TSP), given a set of points in the plane, one seeks 
a shortest tour (closed curve) that visits each point. In the path variant, one seeks a shortest path 
(open curve) that visits each point. If now each point is replaced by a (possibly disconnected) 
region, one obtains the so-called TSP with neighborhoods (TSPN), first studied by Arkin and 
■^J - ! Hassin pQ. A tour for a set of neighborhoods is also referred to as a TSP tour. A path for a set of 

neighborhoods is also referred to as a TSP path. 

For the case of neighborhoods that are (infinite) straight lines, an optimal tour can be computed 
in 0(n 5 ) time [21 El E] (see also [5]), and a \/2-approximation can be computed in 0(n) time [5]. 
For the case of neighborhoods that are (half-infinite) rays, no polynomial time algorithm is known 
for computing an optimal tour, and a \/2-approximation can be computed in 0(n) time [5]. In this 
paper we present linear-time approximation algorithms with improved ratios for these problems. 
The obvious motivation is to provide faster and conceptually simpler algorithmic solutions. As 
mentioned above, while for the case of rays no polynomial time algorithm is known, for the case of 
lines, the known algorithms reduce the problem of computing an optimal tour of the lines to that 
of computing an optimal watchman tour in a simple polygon for which the existent algorithms are 
quite involved and rather slow for large n [2j [8j [9] . 

In this paper we present four improved linear-time approximation algorithms for TSP, for two 
cases of neighborhoods, that are straight lines, and respectively, straight rays in the plane. We 
consider two variants of the problem: that of computing a shortest tour and that of computing 
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Ratio 


Tour (old ratio) 


Tour (new ratio) 


Path (old ratio) 


Path (new ratio) 


Lines 


V2 = 1.41 . . . 


1.28 


2^2 = 2.82... 


1.42 


Rays 


V2 = 1.41 . . . 


1.28 




2.24 



Table 1: Old and new approximation ratios. No approximation for paths on rays was reported in [5]. 

a shortest path visiting the input set. Our algorithms are all based on solving low-dimensional 
linear programs. Our results are summarized in Table [TJ 

Theorem 1. Given a set of n lines in the plane: (i) A TSP tour that is at most 1.28 times longer 
than the optimal can be computed in 0{n) time, (ii) A TSP path that is at most 1.42 times longer 
than the optimal can be computed in 0(n) time. 

For lines, the previous best approximations obtained in linear time were V2 ps 1.41 and 2^2 ps 
2.82, respectively [5]. 

Theorem 2. Given a set of n rays in the plane: (i) A TSP tour that is at most 1.28 times longer 
than the optimal can be computed in 0(n) time, (ii) A TSP path that is at most 2.24 times longer 
than the optimal can be computed in 0[n) time. 

For rays, the previous best approximation for tours was \/2 ~ 1-41 [5] (obtained also in linear 
time, however this was the only approximation known), while for paths there was no approximation 
known. 

Preliminaries. We use the following terms and notations. We denote by x(p) and y(p) the 
x and y-coordinates of a point p. We say that point q dominates point p if x(p) < x(q) and 
yip) < y(q)- For a segment s, Ax(s) and Ay(s) denote the lengths of its horizontal and vertical 
projections. The convex hull of a planar set A is denoted by conv(A). The Euclidean length of a 
curve 7 is denoted by len(7). For a polygon P, let per(P) denote its perimeter. For a rectangle 
Q, let long(Q) denote the length of a longest side of Q. For a ray p, let t{p) denote its supporting 
line. 

The inputs to the two variants of TSP we consider are a set of lines or a set of rays. Let C 
be a given set of n lines, and let T*{£) be an optimal tour (circuit) of the lines in C Let TZ be a 
given set of n rays, and let T*(1Z) be an optimal tour (circuit) of the rays in TZ. 

Following the terminology from (3J [7], a polygon is an intersecting polygon of a set of regions 
in the plane if every region in the set intersects the interior or the boundary of the polygon. The 
problem of computing a minimum-perimeter intersecting polygon (MPIP) for the case when the 
regions are line segments was first considered by Rappaport [7] in 1995. As of now, MPIP (for line 
segments) is not known to be polynomial, nor it is known to be NP-hard. 

Since both lines and rays are infinite (i.e., unbounded regions) finding optimal tours T*{£) 
and T* (TZ) are equivalent to finding minimum-perimeter intersecting polygons (MPIPs) for C and 
TZ respectively. We can assume without loss of generality that not all lines in C are concurrent 
at a common point (this can be easily checked in linear time), thus per(T*(£)) > 0. The same 
assumption can be made for the rays in TZ, thus per(T* (1Z)) > 0. 

The following two facts are easy to prove; see also [H HI \7\. 

Observation 1. If Pi is an intersecting polygon of C, and P\ is contained in another polygon P2, 
then P2 is also an intersecting polygon of C The same statement holds for 1Z. 
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Observation 2. T*(£) is a convex polygon with at most n vertices. Similarly, T*(1Z) is a convex 
polygon with at most n vertices. 

A key fact in the analysis of the approximation algorithms for computing tours is the following 
lemma. This inequality is implicit in [10] ; a more direct proof can be found in [3]. 

Lemma 1. [3l[T0]. Let P be a convex polygon. Then the minimum-perimeter rectangle Q con- 
taining P satisfies per(Q) < - per(P). 

2 TSP for lines 

In this section we prove Theorem [TJ 

TSP tours. We present a ^(l+e)-approximation algorithm for computing a minimum-perimeter 
intersecting polygon of a set C of n lines, running in 0(n) time. If we set e = 1/200, we get the 
approximation ratio 1.28. For technical reasons (see below) we choose e G [1/300, 1/200] uniformly 
at random, and the approximation ratio remains 1.28. The algorithm combines ideas from OHKS]. 
As in [3j, we first use the fact (guaranteed by Lemma[l| that every convex polygon P is contained in 
some rectangle Q = Q{P) that satisfies per(Q) < ^ per(P). In particular, this holds for the optimal 
tours T*{C) and T*(JZ). Then we use linear programming to compute a (1 + ^-approximation for 
the minimum-perimeter intersecting rectangle of C (as in [3]; see also [5]). 

Algorithm Al. 

Let m = ■ For each direction on = i • 2e, i = 0, 1, . . . , m — 1, compute a minimum- 
perimeter intersecting rectangle Qi of C with orientation «j. Return the rectangle with the 
minimum perimeter over all m directions. 

We now show how to compute the rectangle Qi by linear programming. By a suitable rotation 
(by angle a,) of the set C of lines in each iteration i > 1, we can assume that the rectangle Qi 
is axis-parallel. This can be obtained in O(n) time (per iteration). Let {qi,q2,Q3,Qi} be the four 
vertices of Qi in counterclockwise order, starting with the lower leftmost corner as in Figure EJ 
As in [5], let C = C~ U C + be the partition of C into lines of negative slope and lines of positive 
slope. By setting e G [1/300, 1/200] uniformly at random, in each iteration i, with probability 1 
there are no vertical lines in (the rotated set) C. 

Observe (as in [5]), that a line in i G C + intersects Qi if and only if c/2 and are separated 
by I (points on t belong to both sides of £). Similarly, a line in I G C~ intersects Qi if and 
only if q\ and q% are separated by i. The objective of minimum perimeter is naturally expressed 
as a linear function. The resulting linear program has 4 variables x\,X2,yi,y2 for the rectangle 
Qi = [x\,X2] x [yi,2/2] ; and 2n + 2 constraints. 



minimize 2(x2 — x\) + 2(7/2 
y 2 > ax\ + b, 



subject to 



yi < 0x2 + b, 
yi < ax\ + b, 
y 2 > 0x2 + b, 
Xl < x 2 

yi < V2 



Vi) (LP1) 

y = ax + b G C + 
y = ax + b G C + 
y = ax + b G £~ 
y = ax + b G £~ 
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Let Q* be a minimum-perimeter intersecting rectangle of C. To account for the error made by 
discretization, we need the following easy fact; see [H Lemma 2]. 

Lemma 2. |3j. There exists an i E {0, 1, . . . , m — 1} such that per(Qj) < (1 + e) per(Q*). 

By Observations [1] and [21 and by Lemmas [1] and [21 the algorithm Al computes a tour that is 
at most 1.28 longer than the optimal. The algorithm solves a constant number of 4-dimensional 
linear programs, each in O(n) time [6]. The overall time is 0(n). 

TSP paths. The key to the improvement is offered by the following. 

Observation 3. Let Q be a rectangle. Then Q intersects a set of lines C if and only if any three 
sides of Q intersect C 

Proof. Fix any three sides of Q: {si, S2, S3} (each Si is a closed segment). Now if £ is a line 
intersecting Q, then i intersects at least two sides of Q, hence it intersects at least one element of 
{si, S2, S3}, as required. □ 

The next lemma gives a quantitative upper bound on the total length of three shorter sides of 
a rectangle enclosing a curve. 

Lemma 3. Any open curve of length L can be included in a rectangle Q, so that per(Q) — long(<5) < 
\plL. This inequality is the best possible. 

Proof. Let 7 be an open curve of length L = len(7), and let a, b S 7 be the two endpoints of 7. We 
can assume w.l.o.g. that ab is a horizontal segment, and let Q = Q{^) be a minimal axis-aligned 
rectangle containing 7. Write z = \ab\. Let w and h the lengths of the horizontal and vertical 
sides of Q, respectively (i.e., the width and height of Q). It suffices to show that w + 2h < \[2L. 

By construction 7 meets each side of Q. We trace 7 from a to b and subdivide it into a finite 
number of open sub-curves 7^; the endpoints of each sub-curve ji belong to two distinct sides of 
Q. We denote by Si the segment connecting the two endpoints of j{. By concatenating these 
segments (in the same traversal order) we get a polygonal curve connecting a and b. We call this 
(not necessarily unique) curve, a polygonal curve induced by 7; see Figured! 



a b 




Figure 1: Polygonal curves induced by 7 (two examples, in bold). The corresponding rectangles Q are 
drawn in dashed lines. 

For any segment s, we have Ai(s) + Ay(s) < |s|v2: indeed, Ax(s) + Ay(s) = |s|(sina + 
coscv) < |s|\/2, by a well-known trigonometric inequality (here a G [0,7r/2]). By adding the above 
inequalities for all segments Sj (and sub-curves 7^) we obtain 

(Ax(si) + Ay(si)) < (J^ |s f |) V2 < len( 7 )V2 = Ly/2. (1) 
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On the other hand we have ^ Ax(si) > w. indeed, 7 starts at a and meets the left and right 
sides of Q before returning to b, hence the horizontal projections of the segments s, sum up to at 
least 2w — \ab\ = 2w — z > w. Similarly, we have ^ Ay (si) > 2h: indeed, 7 starts at a and meets 
the top and bottom sides of Q before returning to b. Since a and b have the same y-coordinate, 
the vertical projections of the segments Sj cover twice the height of Q. Consequently, we have 

^(Ax( Si ) + Ay( Si )) >w + 2h. (2) 

Putting ([T]) and ((2j) together yields 

w + 2h<J2(Ax( Si ) + Ay(si)) < LV2, (3) 

as required. 

To see that this inequality is tight, let 7 be a two-segment polygonal path made from the two 
unit sides of an isosceles right triangle. Then L = len(7) = 2, while the rectangle Q enclosing 7 has 
sides \/2 and \f2j2 respectively. The lengths of its three shorter sides sum up to ^+2^/2/2 = 
2^/2 = Ly/2. It can be verified that the sum of the three smallest sides of any other enclosing 
rectangle is larger (details in the Appendix), hence the rectangle Q constructed in our proof is 
optimal for 7. The proof of Lemma [3] is now complete. □ 

To compute a TSP path for a set of n lines, we use the algorithm A2 we describe next. This 
algorithm is similar to algorithm Al, described earlier. A2 computes a rectangle in each direction 
from a given sequence. The only difference in the linear program is that instead of minimizing the 
perimeter of an intersecting rectangle, 2(x2 — x\) + 2(y2 — yi), it minimizes the sum of the lengths 
of three sides, namely (X2 — x\) + 2(y2 — Vij- The objective function is not symmetric with respect 
to the two coordinates axes, and so the number of directions m, from algorithm Al, is m = [Jj] 
in algorithm A2. Let now Q* be an intersecting rectangle of C with minimum sum of the lengths 
of three sides. Analogous to Lemma [2] we have 

Lemma 4. There exists an i £ {0, 1, . . . , m — 1} such that 

per(Q 4 ) - long(Qi) < (1 + e) (per(Q*) - long(Q*)). 

By Lemma [3] and Lemma H] the approximation ratio is v2(l + e), and we set e = 1/250 (or 
slightly smaller, as before), to obtain the approximation ratio 1.42. This completes the proof of 
Theorem [TJ 

3 TSP for rays 

As noted in [5] : If the lines are replaced by line segments the problem of finding an optimal tour 
becomes NP-hard. Should the lines be replaced by rays, we get a variant of the problem that lies 
somewhere in between the variant for lines and that for line segments, and whose complexity is 
open. In this section we prove Theorem [2l 

TSP tours. The algorithm Al from Section[2]can be adapted to compute a ^(l+e)-approximate 
tour for a set IZoin rays. Let m = ■ The resulting algorithm A3 for computing an approximate 
tour for n given rays computes a minimum-perimeter rectangle intersecting all rays over all m 
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Figure 2: The rectangle Qi, and two rays, one in TZ± and one in IZ^ that intersect it. 

directions. As before, assume that in the ith iteration, the rectangle Qi = {qi, q%, qz, (74} is axis- 
parallel. A ray in 1Z is said to belong to the ith quadrant, i = 1,2, 3, 4, if its head belongs to the 
ith. quadrant when placed with its apex at the origin. Let TZ = TZi U 7^2 U TZ3 U IZ4 be the partition 
of the rays in 1Z (after rotation) as dictated by the four quadrants. See Figure [2J 
Observe that: 

• A ray p G TZ\ intersects Qi if and only if qi and 54 are separated by £(p), and the apex 
(endpoint) of p is dominated by q^. 

• A ray p G IZ2 intersects if and only if q\ and #3 are separated by £(p), and the apex of p 
lies right and below 54. 

• A ray p G 7^.3 intersects if and only if t/ 2 and 94 are separated by £(p), and the apex of p 
dominates q\. 

• A ray p G 7^4 intersects Qi if and only if gi and (73 are separated by £(p), and the apex of p 
lies left and above 52- 

The constraints listed above correct an error in the old \/2-approximation algorithm from [5], 
where it was incorrectly demanded that the apexes of the rays must lie in the rectangle Qi. Indeed, 
this condition is not necessary, and moreover, may prohibit finding an approximate solution with 
the claimed guarantee of \/2. 

Observe that these intersection conditions can be expressed as linear constraints in the four 
variables, xi, X2,yi,y2- The algorithm A3 computes a minimum-perimeter rectangle intersecting 
all rays over all m directions. For each of these directions, the algorithm solves a linear program 
with four variables and 0(n) constraints, as described above. As such, the algorithm takes O(n) 
time |Bj. The approximation ratio is — (1 + e), and we set e = 1/200 (or slightly smaller, as before), 
to obtain the approximation ratio 1.28. 

TSP paths. We need an analogue of Lemma [2] for open curves. This is Lemma [5] below which 
is obviously also of independent interest. 

Lemma 5. Any open curve of length L can be included in a rectangle Q, so that per(Q) < y/bL. 
This inequality is the best possible. 

Proof. Let 7 be an open curve of length L = len(7), and let a, b G 7 be the two endpoints of 7. We 
can assume w.l.o.g. that ab is a horizontal segment, and let Q = Q{^) be a minimal axis-aligned 
rectangle containing 7. Let w and h be the lengths of the horizontal and vertical sides of Q, 
respectively (i.e., the width and height of Q). 



It suffices to show that 2w + 2h < \fhL. Write w = L/X, where A > 1. The case A = 1 
corresponds to a degenerate enclosing rectangle, when 7 is a line segment, and for which 2w + 2h = 
2w = 2L < VbL. We can therefore assume in the following that A > 1. By construction 7 meets 
each side of Q. Arbitrarily select a point of 7 on each of these sides to obtain a polygonal open 
curve 71 connecting a and b and passing through these intermediate points (and still enclosed in 
Q). By construction, the intermediate points are visited in the same order by 7 and 71. By the 
triangle inequality, len(7i) < len(7). 

Successively reflect the rectangle Q with respect to the sides containing the intermediate points 
in the order of traversal. See Figure [3] for an illustration. Let b' be the final reflection of the end- 
point b of 7. The segments composing 71 can be retrieved in the reflected rectangles; they make 
up a polygonal path 72 connecting a and b' . By construction we have len(7i) = len(72). 




b 



Figure 3: The reflection argument (two examples, in bold). The small shaded triangles indicate the lower 
left corner of Q in all subsequent reflections. 

It is easily seen that Ax(ab') = 2w — \ab\ and Ax(ab') = 2h. By Pythagoras' Theorem, 

Ax 2 (ab') + Ay 2 {ab') = \ab'\ 2 . 
Since the shortest distance between two points is a straight line, len(72) > \ab'\. It follows that 

L 2 = len 2 ( 7 ) > len 2 ( 7 i) = len 2 ( 72 ) > |a6'| 2 = (2w - \ab\) 2 + Ah 2 . (A) 
Obviously, \ab\ < w, thus from @ we deduce that 

w 2 + Ah 2 < L 2 . (5) 
Substituting w = L/X in §5§ yields L 2 /A 2 + Ah 2 < L 2 , hence 



It follows that 



per(Q) = 2„ + 2h < ( f + ) L = ( ~ ^ ) L. (6) 



Consider the function 
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Its derivative is 



fl(x) = 7^t-2-VA^T = A 2 -2VA^T-(A 2 -1) = 1 - 2y^T 
[ ' A 2 A 2 v / A 2 ^T A 2 v / A 2_r T 

It can be easily verified that /'(A) vanishes at A = \/5/2, and that /(A) is increasing on the 
interval (l,\/5/2] and decreasing on the interval [\/5/2, oo). Hence /(A) attains its maximum at 
A = y/E/2, that is, /(A) < f(VE/2) = y/h. According to ©, we have per(Q) < Ly/b, as desired. 

To see that this inequality is tight, let 7 be a two-segment polygonal path made from the two 
equal sides of an isosceles triangle with sides 1, 1, and 4/\/5. In this case, L = 2 and per(Q) = 
2(4/^5 + l/>/5) = 2y/E = Ly/E. It can be verified that the perimeter of any other enclosing 
rectangle is larger (details in the Appendix), hence the rectangle Q enclosing 7 constructed in our 
proof has minimum perimeter. The proof of Lemma [5] is now complete. □ 

Given TZ, we compute an approximation of the optimal TSP path by using algorithm A3. Let 
now 7 be an optimal TSP path for the rays in TZ, where L = len(7). Since every ray in TZ intersects 
7, every ray in 1Z intersects the rectangle Q = Q(j), as defined in the proof of Lemma For 
each of the m directions, the algorithm A3 computes a minimum-perimeter rectangle (of that 
orientation) intersecting each ray in 7Z. Thus A3 computes a rectangle (i.e., a closed path) that 
is a (1 + e)-approximation of the minimum-perimeter rectangle intersecting each ray in 1Z. For 
e = 1/1000, its perimeter is at most (1 + e)\/EL < 2.24L, as claimed. This completes the proof of 
Theorem [2j 

4 Final remarks 

Interesting questions remain open regarding the structure of optimal TSP tours for lines and rays, 
and the degree of approximation achievable for these problems. For instance, the following two 
problems (one new and one old) deserve attention. 

(1) Is there a polynomial-time algorithm for computing a shortest tour (or path) for a given set 
of rays in the plane? 

(2) Can one compute a good approximation of a shortest TSP tour (or path) for a given set 
of lines in 3-space? Note that this variant with parallel lines reduces to the Euclidean TSP 
for points in the plane (namely the points of intersection between the given lines and an 
orthogonal plane), so it is NP-hard. See also [4]. 
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Appendix 

A tight bound for Lemma [3], Let Q be a minimal rectangle containing 7 = acb, whose width 
and height are w and h, respectively. Let A = conv(7) = conv(ac6), where Zacb = tt/2. It suffices 
to show that w + 2h > L\[2L = 2\/2- By the minimality of Q, at least one vertex of A must 
coincide with a corner of Q, say the lower left corner q\. If c = q\, then Q is a unit square, thus 
w + 2h = 3 > 2\/2- Assume now that a = q±, as in Figure [U 



We have w = sin a + cos a and h = cos a, where a G [0, 7r/4]. Hence w + 2h = (sin a + cos a) + 
2 cos a = 3 cos a + sin a. Consider the function f(a) = 3 cos a + sin a, where a G [0, 7r/4]. Its 
derivative, f'(a) = —3 sin a + cos a vanishes at a = arctan(l/3), and is positive on (0, arctan(l/3)) 
and negative on (arctan(l/3), vr/4). Hence f(a) attains its minimum at one of the endpoints of 
the interval [0,vr/4]. We have /(0) = 3 and /(tt/4) = 4^/2/2 = 2^2, therefore f(a) > 2^2 for 
a G [0,71-/4]. 

A tight bound for Lemma [5j Again, let Q be a minimal rectangle containing 7 = acb, whose 
width and height are w and h, respectively. Let A = conv(7) = conv(acft), where Zacb = 2 arctan 2. 
It suffices to show that 2{w + h) > L\fhL = 2yE. By the minimality of Q, at least one of the 



c 




a 



Figure 4: A rectangle containing 7. 



9 




two vertices a and b of A must coincide with a corner of Q, say the lower left corner q\. We have 
a = qi, as in Figure [5j We distinguish two cases: 

Case 1. The vertex b lies on the right side of Q, and c lies on the top side of Q, as in 
Figure [SJdeft). We have w = -^gcosa and h = -^sina + sin/3, where a E [0, arctan(l/2)], and 
/3 = 7r/2 — (-7T — (7r/2 — a) — arctan(l/2)) = arctan(l/2) — a. This yields 

sin /3 = sin(arctan(l/2) — a) = sin arctan(l/2) cos a — cos arctan(l/2) sin a 

1 2 
= — = cos a = sin a. 

It follows that 

/ 4 4 1 2 . \ / 2\ 

2(u; + ft ) = 2 cos a H = sin a H = cos a = sin a =2 v 5 cos a H = sin a . 

V5 VE J V VE J 

Consider the function /(a) = \/5cosa -f (2/\/5)sina, where a E [0, arctan(l/2)]. Its derivative, 
/'(a) = — \/5sinQ + (2/-v/5)cosa vanishes at a = arctan(2/5), and is positive on (0, arctan(2/5)) 
and negative on (arctan(2/5), arctan(l/2)). Hence f(a) attains its minimum at one of the end- 
points of the interval [0, arctan(l/2)]. We have /(0) = V5 and /(arctan(l/2)) = 2 + 2/5 = 12/5 > 
V§, therefore f(a) > y/b for a E [0, arctan(l/2)]. Equivalently, per(Q) = 2(w + h) > 2\/5 in 
this case. (It may be noted that the rectangle with perimeter 2 • 12/5 = 24/5 corresponding to 
a = arctan(l/2) is flush with the unit side cb.) 

Case 2. The vertex b coincides with the upper right corner of Q, and c lies in the interior of 
Q, as in Figure [5j right). By symmetry, we can assume that w > h. We have w = (4/v5)cosa) 
and h = (4/^5) sina, where a = (3 + arctan(l/2), and f3 E [0, 7r/4 — arctan(l/2)]. This yields 

4 

w + h = ^(cos(/3 + arctan(l/2)) + sin(/3 + arctan(l/2))). 
V5 

The above expression attains its minimum at f3 = 0, thus 

g 

per(Q) = 2(w + h) > —= (cos arctan(l/2) + sinarctan(l/2)) 
" V5 

— (— + — ^ - — — - — > 2V5 

~ \/5 \V5 V5J~V5 V5 ~ 5 > 

in this second case. 

We therefore always have per(Q) > 2y5, as claimed. 
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